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FOREWORD 


This  work  was  performed  to  provide  a  design  tool  for  use  in  estin  ng  the 
aerodynamics  of  guided  and  unguiJed  projectiles.  Support  for  the  \vcrK  w„  provided 
by  the  Naval  Ordnance  Systems  Command  under  ORDTASK  35A-53I  «0-l/UF 
32-323-505. 

The  current  report  was  reviewed  and  approved  by  Mr.  D.  A.  Jones,  ill.  Head 
of  the  Aeroballisties  Group  and  by  Mr.  C.  A  Cooper,  Head  of  the  Guid  .1  Projectile 
Division. 
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ABSTRACT 

Several  theoretical  and  empirical  procedures  are  combined  into  a  Single 
computer  program  to  predict  lift  drag,  and  center  of  pressure  on  .|uite  general 
wing-body  geometries.  The  method  is  applicable  for  Mach  number  zero  to  three  and 
angle-oi-attack  zero  to  about  fifteen  degrees.  Computed  results  for  several 
config  t  ratio  ns  compare  well  with  experimental  and  other  analytical  results.  It  costs 
about  five  dollars  per  Mach  number  to  compute  the  static  aerodynamics  of  a  typical 
wing-body  shape  on  the  CDC  6700  computer. 
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I.  INTRODUCTION 


The  goal  of  the  present  research  is  to  develop  the  capability  to  compute  static 
aerodynamics  on  configurations  such  as  guided  and  uuguided  projectiles  for  the  Mach 
number  range  zero  to  three  and  angle-of-attaek  range  zero  to  about  twenty  degrees. 
The  Mach  number  and  angle-of-attack  range  cover  a  majority  of  present  and 
probable  future  design  requirements  for  gun  launched  weapons.  This  work,  therefore, 
is  a  natural  extension  to  the  body  alone  aerodynamic  prediction  methodology 
developed  in  Reference  1.  Included  herein  is  an  outline  of  the  theoretical  and 
empirical  methods  used  to  compute  lift.  drag,  and  pitching  moment  on  wings  and 
tails.  Also  included  arc  procedures  to  determine  the  various  interference  effects 
which  occur  between  the  wing-tail  or  body-tail.  A  description  and  listing  of  the 
entire  wing-body  computer  program  will  appear  as  part  II  of  this  report. 

The  overall  guiding  principle  of  the  current  work  is  to  use  anah  tical  methods 
which  yield  reasonable  accuracies  and  icqutre  reasonable  computational  time.  For 
areas  where  the  state-of-the-art  is  such  that  analytical  procedures  do  not  meet  the 
above  requirements  on  accuracy  and  cost,  empirical  procedures  are  employed.  This 
problem  occurs  mainly  in  transonic  flow. 

There  have  been  several  works  which  pertain  to  the  present  problem.  The  first 
of  these  is  that  of  Saffell,  ct  ai‘ ’ 1  who  devoir,  ed  a  computer  program  to  compute 
static  aerodynamics  on  low  aspect  ratio  missile  configurations.  The  method  could  lu 
applied  for  large  angles-of-attack  and  for  subsonic  thiough  supersonic  Mach  numbers. 
However,  the  drag  was  computed  using  handbook  techniques  and  the  lift  of  the 
wing  alone  was  found  from  an  empirical  formula  for  low  aspect  ratio  wings.  Thus, 
drag  results  aie  quite  inaccurate  tat  small  angles-of-attack),  as  well  as  lift  for  high 
aspect  ratio  configurations. 

Another  method  which  is  available  for  spin-stabilized  (or  body  alone)  projectiles 
is  the  GK  'Spinner'  program.* 3  *  This  program,  which  is  completely  empirical,  gives 
good  results  for  most  stand;. *d  shaped  projectiles  but  is  too  limited  in  scope  to 
meet  the  present  needs. 

The  most  detailed  method  of  those  previously  available  is  that  of  Woodward.14 1 
Woodward  uses  perturbation  theoiy  to  compute  the  pressure  distribution  on 
wing-body  combinations  in  subsonic  and  supersonic  flow.  However,  the  bodies  must 
be  pointed  and  the  wing  leading  edge  sharp  Also,  he  does  not  calculate  the  base 
and  skin-friction  drag  or  the  non-linear  angle-of-atiack  effects. 

As  is  apparent  irom  the  above  discussion,  none  of  the  piewmis  works  can  in 
itself  accurately  compute  total  lift.  drag,  and  pitching  moment  on  wmg-body 

1 


combinations  for  Mach  numbers  0  to  3,  and  angles-of-attack  0°  to  20°.  Moreover, 
no  attempt  is  made  to  handle  the  complicated  body  and  wing  geometries  (see 
Figure  1)  which  arise  fiom  gun  launched  guided  weapons. 
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a.  SPIN  STA3IUZED  PROJECTILE 


b.  UNGUIDEO  RN-  STABILIZED  PROJECTILE 


FIGURE  l  BASIC  CONFIGURATIONS  ( ILLUSTRATION  ONLY! 


U.  ANALYSIS 


The  approach  of  the  piesent  work  is  to  break  the  missile  configuration  down 
into  its  individual  components  composed  ol  body  alone,  wines,  and  canards  and 
then  to  account  separate!;,  for  the  various  interference  effects.  1  Ins  is  opposed  to 
mathematically  modelling  the  entire  configuration  simultaneously  as  w.»s  done  in  the 
work  of  WooUwa»d<4>  Woodward  considered  the  wing-body  simultaneously  through 
an  appropriate  source  and  sink  distribution,  where  the  sources  from  the  body  were 
allowed  to  intluence  the  wing  solution  and  vice  versa.  This  is  (airly  straightforward 
for  pointed  nosed  bodies  .0  d  wings  with  sharp  leading  edges,  but  for  blunt  nosed 
bodies  and  blunt  leading  edge  wings,  as  is  the  case  in  the  present  work,  the 
problem  is  complicated  considerably.  Moreover,  assuming  perturbation  theories  are 
used  to  calculate  the  imiscid  aerodynamics,  the  approach  ol  considering  the 
individual  components  separately  should  yield  total  forces  and  moments  which  aie  as 
good  as  those  that  would  be  obtained  by  considering  the  configuration  as  a  whole. 

A.  Body  Alone  Aerodynamics 

The  body  alone  aerodynamic  analysis  appears  in  Reference  f  and  will  not 
be  repeated  here.  However,  a  summary  of  the  various  methods  for  computing  body 
alone  aerodynamics  appears  in  Figure  2.  All  the  methods  ate  slant  rd  in  tin* 
literature  (References  5  through  10)  with  the  exception  of  the  empirical  schemes 
derived  for  transonic  lift  and  wave  drag  and  the  combined  Nevvtonian-perturba.ion 
theory  for  calculating  nose  wave  diag  in  supersonic  How  The  combined 
Newtonian-, H’rturbation  theory  was  developed  so  icasomible  icsults  for  static 
aerodynamics  could  be  obtained  at  low  supersonic  Mach  mmitvis  tor  blunt  nosed 
configurations.  Previously,  the  available  theories  were  eithei  too  mac  curate' 1 1 1  or  too 
complicated"’’  for  use  at  the  lower  end  of  the  supersonic  Mach  range 
(1.2  <  M.  <  2.0). 

B.  Tail  or  Canard  Drag 

As  m  the  ease  of  an  :.\is\  mmcirie  body,  tile  three-dimensional  wing  ding 
is  composed  of  wave,  skin-friction,  and  a  trailing  edge  separa'ion  drag  if  the  trading 
edge  is  blunt  or  the  rear  section  of  the  an  foil  has  a  large  slope  lius  trailing  edge- 
separation  drag  is  analogous  to  base  pressure  drag  on  a  body  <>t  (evolution.  In 
addition  to  these  drag  components  the  wings  can  cause  an  additional  drag  from 
body  base  pressure  changes  due  to  the  presence  of  fins.  P'aeh  ol  these  wing  drag 
components  will  be  treated  .cpaiatelv  below 
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METHODS  USED  TO  COMPUTE  BODY  ALONE  AERODYNAMICS 


1.  Wave  Drag 


It  will  be  assumed  apriori  that  the  wing  is  symmetric  at  nit  the  x-y 
plane  so  that  no  „„.T.ber  is  present.  Furthermore,  the  wing  will  be  assumed  to  be 
th  n  wnh  cither  a  modified  double  wedge  (Figure  3A)  cr  biconvex  airfoil  section 
(Figure  3B).  However,  by  assuming  a  modified  double  wedge  only  requires  that 
straight  lines  exist  between  points  A  and  B,  B  and  D,  and  D  and  E  (see 
Figure  3A).  These  straight  lines  could  then  be  any  percentage  of  the  entire  chord. 
For  example,  if  BD  were  zero  the  airfoil  would  be  a  double  wedge  design  or  if 
both  BD  and  DE  were  ze.o,  the  airfoil  would  be  a  wedge.  Also,  either  the  biconvex 
oi  modified  double  wedge  design  may  have  blunt  leading  and  trailing  edges  and  the 
thickness  to  chord  ratio  may  vary  along  the  span.  The  wing  ge  lerators  GK,  HL,  IM, 
and  JN  are  assumed  straight. 

Since  the  wing  is  thin,  the  linearized  three-dimensional  equation  of 
motion  governing  the  steady  flow  field  is:(13) 


"  4>yy  "  =  0 


(1) 


where  subscripts  indicate  partial  differentiation.  Here,  the  velocity  potential  0  is 
i  da  ted  to  the  perturbation  velocities  by: 


0X  =  u 


(2a) 


0y  =  V  (2b) 

0,  -  w  (2c) 


The  boundary  conditions  required  for  the  solution  of  the  linear  partial  differential 
I  quation  (!)  are  that  the  fi  -  must  be  tangent  to  the  surface: 


w(x.y)  =  0  (x.y.O  )  =  0  (x,y,0+)  =  — (x,y)  (3) 

1  Ox 

and  that  the  perturbation  velocities  must  vanish  upstream  of  the  most  forward  point 
ol  the  wing.  Referring  to  Figure  3A,  this  most  forward  point  is  at  x  =  O'  so  that 
the  second  boundary  condition  is: 
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(4) 


N*  * 

t 


-  C\  10  =  d>  (O.y.z)  =  0 


Hie  plus  or  minus  superscript  means  the  pjirticula.  axis  is  approached  Irom  the 
positive  or  negative  side,  respectively  Due  to  the  symmetrv  of  the  airloii. 
kquation  (3i  indicates  that  it  makes  no  difference  irom  which  side  one  appioaches 
the  axis  /  =  0. 

liquation  ( 1 )  is  valid  only  where  the  perturbation  velocities  aie  small. 
This  means  that  in  the  neighborhood  of  a  blunt  leading  or  trailing  edge,  some  other 
method  must  be  applied.  Consider  first  a  gem  kiI  three-dimensional  wing  with  sharp 
leading  and  trailing  edges. 

The  general  solution  to  hquation  (I)  along  the  airfoil  surface  (/  £  0) 

is"1' 


0(\.y  0) 


w(x,.y,  )dx1dy1 
yTx-x,  )--p’>( y-  y,  )2 


where  Z  indicates  the  region  of  integration.  The  source  strength  \v(x, ,  y,)  is  related 
to  the  'oeal  slope  of  the  airfoil  surface  through  the  boundary  vendition 
liquation  (3). 


In  previous  works,  w( \ ,  )  was  assumed  constant  or  a  Inaction  of  \ 
only  (the  slope  of  the  airfoil  surface  was  the  same  all  along  the  span).  so  the 
integration  of  the  above  integral  could  be  carried  out  in  closed  form  for  simple 

wing  planforms.' 1  !  1  In  the  present  analysis,  the  slope  of  the  wing  is  allowed  i<» 

vary  in  the  spanwise  as  well  as  the  chord  wise  direction  so  the  in  egialiou  oi 

I  qualion  (5)  cannot,  in  general,  he  carried  out  in  closed  form.  The  most 

straightforward  method  of  solution,  is  then  to  define  the  slope  of  the  given  surface 
and  carr>  out  the  double  integration  In  numerical  quadiature.  However,  one  must 
be  aware  ot  the  singular  nature  of  the  double  integi.il  where  (x~  x(  )’  =  p’t\ 
during  the  integration  piocess.  Anothei  alternative  is  tc»  assume  that  on  a  small 
element  of  the  wing  surface,  \vt\..\,)  is  constant.  Iihh.  if  the  region  ot  tukgintK  n 
of  lquation  (5).  Z  is  assumed  to  be  ovet  a  small  element  of  the  wing,  one  nui.v 
write. 
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Equation  (/>»  is  now  in  the  form  given  m  Reference  13  for  simple  planfonn 

geometries  and  the  integration  can  be  carried  out  directly.  Again  it  should  be 

emphasized  that  w(x.y)  is  the  slope  of  the  airfoil  surface  at  a  given  point  and  vanes 
for  each  element  on  the  wing. 

The  closed  form  solution  of  Equation  (6)  depends  on  whether  the 
wing  generators  are  subsonic  or  supersonic.  Referring  to  Figure  3A,  a  subsonic  wing 
generator  waul.!  exist  if  the  Mach  number  normal  to  line  GK  were  less  than  oil* 
and  a  supersonic  wing  generator  would  exist  if  Mn  >  1.  The  same  applies  to  each 
of  the  other  wing  generators  HL,  IM,  and  IN.  For  the  biconvex  airfoil.  Figure  313,  a 
continuous  distribution  of  wing  generators  is  placed  between  the  lines  GK  and  JN. 
Each  wing  generator  is  analogous  tc  a  line  of  sources  and  sinks  with  strengths 

-sufficient  to  keep  the  flow  tangent  to  the  surface. 

a.  Subsonic  Source  or  Sink  Line  (SOSL) 


If  the  wing  generator  is  subsonic,  the  induced  velocity  at  a  given 
point  P.  due  to  the  SOSL,  is  dependent  on  the  location  of  P  relative  to  the  SOSL. 
Referring  to  Figure  4A.  if  P  -  P,  the  induced  velocity  is:l,,) 


0 


—H^^cosh-' 

rr/Vr?-  '  1 


(7) 


whc'V  w  is  determined  from  the  boundary  condition  and  is  (for  the  airfoil  section 

Jl  y  =  w 


dz 


MM 


In  Equation  (7),  the  definitions 


V  ~  K//S 


k  =  tan  A 


a  = 


kv\ 


1 7a  I 


1;  ive  been  used  If  P 


IV  the  induced  velocity  at  P,  due  to  a  given  SOSL  is: 
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FIGURE  4A.  TRIANGULAR  SURFACE  SYMMETRIC  ABOUT  X  AXIS 


I  WING  TIP 


FIGURE  4B.  WING  TIP  EFFECTS 


FIGURE  4.  SUBSONIC  SOURCE  OR  SINK  LINE 


2wtxl,*v..,)  ,  .  hr  o'- 

0,  ~  - - cosh  / - " 

'  -  1  v  I  -  <r 


\i  the  wing  tip.  there  is  an  additional  disturbance  within  the  Mach  line  emanating 
•Yom  the  tip  leading  edge  (Figure  4B).  The  induced  velocity  in  this  region.  P  -  P> 


_  v/’xer-Vfj2  .  n~  +  iff! 

'  nfis/rr'-  1  "* S'  t?(!(»1+  1) 


The  absolute  value  of  a  is  taken  because  o  is  actually  negative  for  the  point  P.. 
I  he  induced  velocity  at  any  point,  say  P  =  .  outside  of  the  Mach  lines  emanatine 
Tom  the  beginning  ol  the  SOSl.  is  zero  since  this  point  is  out  of  the  /one  i*f 
influence. 

b.  Supersi  nie  SOSL. 

If  the  wing  generator  is  supersonic,  the  Mach  lines  from  point  0 
in  Figure  5A  lie  behind  the  SOSL.  If  in  Figure  5.  P  r-  P, .  then  the  unliked 
velocity  at  P,  due  to  the  disturbance  caused  by  the  SOSL  is. 1  ^ 


w(\„  s  ,yp ,  > 
tV  I  ~ 


If  P  =  P,.  the  induced  velocity  is 


w(W 

JTp/r-  a- 


7r  -  2  sin  1 


I  ii2  -  j2 


Keferring  to  Figure  5B,  the  additional  induced  velocity  inside  the  area  bounded  by 
the  tip  and  the  Mach  line  emanating  from  the  tip  (P  =  P3)  is: 


wtv  ,y.  )  f  , 
r3  p3  , lie; 

~~  j-j—— tOS  . 

7r/5v  1  if  [_i?tl+|on 
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FIGURE  58.  WING  TIP  EFFECTS 


FIGURE  5.  SUPERSONIC  SOURCE  OR  SINK  LINE 
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Again  if  P  =  P4.  the  point  is  out  of  the  zone  of  influence  of  the  SOSL  and  thus 
the  induced  velocity  is  zero. 


The  induced  velocity  at  a  given  point  on  any  wing  geometry  can 
now  be  computed  by  the  proper  superposition  of  the  triangular  SOSL  shown  in 
Figures  4  and  5.  This  is  because  of  the  linear  nature  of  the  governing  flow-field 
liquation  (1).  As  an  example  of  how  the  above  superposition  principle  works, 
consider  the  wing  shown  in  Figure  6.  For  simplicity,  the  slopes  Xj  and  X2  arc 
constant.  The  ving  AHJD  can  be  represented  by  the  superposition  of  five  SOSL. 
The  first  has  the  plan  form  AEH  and  source  intensity: 


W<VV  =  V.x, 


where  xt  <s  the  slope  of  the  segment  AB.  The  second  has  tiie  planfonn  BJF  and 
intensity 


W(VV  =  (x?. "  X1  )V« 


and  the  third  the  plan  form  DJG  and  intensity 


w(xp.yp)  =  -x2Vw 


The  other  two  SOSL  represent  the  tip  effects.  They  are  the  planforms  11JL  and  UK 
and  have  source  intensities  of  opposite  signs  than  those  representing  the  wing. 

The  above  procedure  can  be  applied  to  a  wing  of  general 
plan  form.  The  only  difference  is  that  for  each  point  in  q«  estion,  the  slope  is  not 
constant  as  was  the  ea>e  in  the  simplified  example.  Then  for  some  general  point 
located  on  the  wing  surface,  the  total  induced  velocity  due  to  all  sources  and  sinks 
is  found  by  applying  one  of  the  Equations  (7)  through  (12)  for  each  SOSL.  The 
particular  equation  applied  depends  upon  the  location  of  the  point  relative  to  tile 
SOSL  and  the  Mach  line  as  discussed  earlier.  These  individual  contributions  are  then 
summed  to  get  the  total  induced  velocity.  Knowing  the  total  induced  velocity  at  a 
point  allows  one  to  calculate  the  pressure  coefficient  at  the  given  point  by: 
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,V 


Cp!x.y)  =  -2ov(x,y,0) 


(13) 


The  pressure  coefficient  can  be  calculated  at  a  given  number  of  spunwise  and 
ehordwise  locations.  The  drae  of  a  given  airloil  sc  'lion  at  the  spanwise  station 
y  =-■  yA  is  then 


(14) 


The  total  drag  for  one  'in  of  semispan  b/2  is  then: 


C, 


j  rb  i  2 

5  J  C* 

w  4. 


c(y)dy 


(15) 


where  S„ 


=  b/2(e  +e().  For  cruciform  fins,  the  total  drag  coefficient  is: 


rb  2 

r'>  ■  t .  r 


:dc(y)dy 


( 1(>) 


If  it  is  desired  to  base  the  drag  coefficient  on  the  body  cross-sectional  area,  the 
Foliation  (16)  most  be  multiplied  by  (he  factor  Sw/S  .. 

liquations  (14)  and  (!(>)  can  be  integrated  by  numerical 
quadrature  it  the  generators  of  the  wing  surface  are  supersonic.  If  the  generators  ate 
subsonic,  linear  theory  indicates  the  pressure  coefficients  go  to  infinity  at  the  wing 
generators.  Physically,  this  cannot  be  true  which  means  that  for  subsonic  SOSL,, 
linear  theory  is  not  valid  at  the  SOSL.  The  reason  is  that  the  velocity  perturbations 
in  the  vicinity  of  the  discontinuities  are  no  longer  small,  violating  one  of  the 
assumptions  in  linear  theory.  However,  the  velocity  perturbations  are  small  a  slight 
distance  from  the  SOSL  so  that  linear  theory  can  be  applied.  Numerical  experiments 
Indicated  a  distance  ol  liv;  thousandths  of  the  chord  length  from  the  SOSL  was 
sufficient  and  the  value  of  pressure  calculated  at  this  point  was  assumed  to  exist  up 
to  the  SOSL. 


The  previous  analyst  applies  to  airfoils  with  sharp  leading  and 
•railing  edges.  If  the  airfoil  leading  or  trailing  edge  is  blunt,  some  other  method 
must  be  applied  in  the  vicinity  of  the  blunted  portion  because  the  assumptions  ol 


It 


perturbation  theory  are  violated  there.  In  analogy  to  the  work  of  Reference  1, 
modified  Newtonian  Theory  will  be  applied  to  the  blunt  leading  edges  and  an 
empirical  afterbody  separation  pressure  correction  applied  at  the  blunt  trailing  edges  is  will 
b,  discussed  later. 


The  modified  Newtonian  pressure  coefficient  is 


Cp 


(17) 


where  0  is  the  angle  between  a  tangent  to  the  local  body  surface  and  the 
freest  ream  direction  and  where  the  stagnation  pressure  coefficient  behind  a  normal 
shock  is: 


•yM; 


* 
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r  •>+! 

yn 

(7  +  1  )M“ 
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> 

bvMi  -  <7~  1) 

1 

t 

(18) 


If  the  blunt  leading  edge  of  the  wing  is  cylindrical  in  a  direction  perpendicular  to 
the  leading  edge,  then  this  circular  shape  appears  as  an  ellipse  in  the  streamwi.se 
direction  for  sweptbaek  wings.  Thus,  for  a  given  point  on  the  airfoil  leading  edge 
with  coordinates  (x.y.z),  it  can  be  shown  that  (see  Appendix  A): 


0(:\y.z.)  =  tan* 1 


cos  A  i 


rt  ,  (y)~  x  cos  A, 


(19) 


Note  that  f-quation  (llh  assumes  the  leading  edge  radius  may  vary  along  the 
span,  that  is  rf  ,,  =  r,  ,.(>).  Ihe  pressure  coefficient  over  the  elliptical  leading  edge 
can  now  he  calculated  at  each  airfoil  section  by  combining  liquations  (17),  (18), 
and  (19). 


The  question  that  remains  is  where  does  one  start  the 
perturbation  theory  aft  of  the  btunt  leading  edge  and  where  does  one  end  the 
modified  Newtonian  Theory.  Before  proceeding  to  answer  this  question,  it  is  helpful 
to  review  the  work  of  Reference  1,  which  combined  the  second  order  perturbation 
theory  of  Van  Dyke  with  the  modified  Newtonian  Theory  to  calculate  wave  drag  on 
blunt  bodies  ot  revolution.  In  that  work,  the  perturbation  theory  was  started  as  far 
Upstream  on  the  spherical  cap  as  possible  while  still  getting  reasonably  accurate 
perturbation  pressure  coefficients.  It  was  found  that  slopes  of  25-30  degrees  were 


optimum.  Although  the  results  were  unpublished,  it  was  found  in  that  woiJc  that  if 

first-order  perturbation  theory  were  used,  this  angle  must  be  reduced  'to  about  15 

degrees.  Also,  it  was  found  that  second -order  theory  accounted  very  wcil  1  ir  the 
over-expansion  region  around  the  spherical  cap  whereas  the  first-order  theory  dfc. 
not.  The  important  analogy  to  be  drawn  from  this  discussion  is  that  for  th.ee 
dimensional  wings,  a  first-order  (rather  than  a  second -order)  theory  is  combined 
with  modified  Newtonian  Theory  to  calculate  wave  drag  when  the  leading  edge  is 
blunt.  Thus,  in  analogy  to  bodies  of  revolution,  one  would  intuitively  expect  the 
angle  where  perturbation  theory  begins  to  be  around  15°.  A  discontinuity  in 
pressure  coefficients  of  Newtonian  and  perturbation  theory  is  expected  at  the  match 
point  due  to  the  failure  of  the  first-order  theory  to  account  for  the  over-expansion 
region. 

Figure  7  compares,  qualitatively,  the  first  and  second -order 
theories  when  combined  with  modified  Newtonian  Theory.  The  second-order  theory 
is  started  at  say  0n  =  30°  and  the  flow  allowed  to  overexpand  around  the  shoulder 
and  recompress  downstream.  The  pressure  coefficient  will  usually  intersect  that  of 
Newtonian  Theory  and  this  is  the  match  point,  x,.  indicated  on  the  figure.  This 

combined  theory  usually  follows  the  experimental  data  reasonably  well.  The 

first-order  theory  starts  at  0{  -  15°  and  the  pressure  coefficient  does  not  usually 
decrease  enough  to  intersect  the  Newtonian  theory.  Hence,  a  discontinuity  in  the 
pressure  coefficient  curve  for  the  entire  configuration  exists  at  the  match  point,  x,. 
For  drag  calculations,  the  Newtonian  Theory  is  used  up  to  the  point  xt  (which 
usually  will  be  at  the  juncture  of  the  cylindrical  leading  edge  and  the  airfoil  after¬ 
body)  and  perturbation  theory  past  x,.  The  drag  coefficient  of  the  blunt  leading 
edge  is  (see  Appendix  B): 


whore 


The  drag  of  the  section  aft  of  the  cylindrical  leading  edge  can  be  found  by 
numerical  quadrature  of  Equations  (14)  and  (1(>). 


(ri.  i  \  +  (rLl-  \ 


4  Rnv,bCl>t>  cos2  A 
Sn-» 


co) 


18 


2.  Skin-Friction  Drag 

In  general,  the  boundary  layer  consists  of  a  laminar,  transitional,  and 
turbulent  region;  however,  no  appreciable  errors  in  total  drag  will  be  realized  by 
assuming  a  laminar  and  turbulent  region  with  transition  taking  place  instantcously. 
Moreover,  transition  from  laminar  to  turbulent  How  is  assumed  to  occur  at  a 
Reynolds  number  of  500,000  based  on  the  wing  mean  aerodynamic  chord. 

According  to  Van  Driest,  the  skin-friction  drag  coefficient  for  a  single 

wing  is  then: 


1 

f 

r  1 

°Af  - 

j(C|r)r  + 

{ 

<Cr„ );  -  (CrT )- 
L  J 

^  4 

Here  (C^,)-  and  (C,T )-  are  the  mean  skin-friction  coefficients  ol  turbulent  flow 
based  on  the  mean  chord  and  transitional  location  on  the  mean  chore.,  respectively. 
They  are  computed  by  solving: 


0.242 


A(CfT)»(Tw/T.  )' 


~  (sin  1  B.  +siir' 


implicitly  tor  CV^  where  t?  for  air  is  0.76.  The  constants  A,  B, ,  and  B,  of 
Equation  (22)  are  defined  by: 


13  =  ^ - : — 

1  ( F2  +  4 A2  )v‘ 

r  n  '4 

(T~  DMj 
2T/T 


B2  ~  (F-  +  4A2  )'4 

F  =  1+  (T  -  1  )/2Mj 

T  /T 

W 1  eo 


The  Reynolds  Number  of  Equation  (22)  is 


V 

Rn  =  « 


(23) 
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where  i  is  either  c  or  x.  The  wall  temperature  can  be  related  to  the  freestream 
temperature  by1: 


r  ,'T  =  1  +  0.9^-  M2 

V.  <  / 


(24) 


Finally,  the  laminar  skin-friction  coefficient  C,  is0 

H 


(C,  )r  = 


1.328 


n/*n 


(25) 


where  the  Reynolds  number  is  based  on  x. 

3.  Trailing  Edge  Separation  Drag 

If  the  trailing  edge  is  blunt  or  if  its  slope  is  large,  the  boundary  layer 
will  separate  somewhere  on  the  rear  of  the  wing.  This  results  in  a  high  drag  region 
similar  to  that  on  the  base  of  a  projectile,  except  here  the  separation  is  a 
two-dimensional  as  opposed  to  a  three-dimensional  phenomenon.  The  pressure  on  the 
rear  of  the  wing  will  then  be  that  of  a  two-dimensional  rearward  facing  step. 
Chapman* 151  presents  experimental  results  for  a  blunt  wing  with  no  slope  at  the 
trailing  edge.  These  results  are  presented  in  Figure  8  as  a  function  of  Mach  number. 
Note  that  the  data  for  Mx  <1.1  has  been  extrapolated  based  on  the  general  shape 
of  the  three-dimensional  base  pressure  curve  presented  in  Reference  1.  If  (rrI  )(  and 
(ij!  ),  are  the  radius  of  the  trailing  edge  hlcntnoss  at  the  root  and  tip.  respectively, 
the  trailing  edge  separation  drag  for  one  fin  is  (see  Appendix  C): 
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I  oi  cruciform  fins  this  become' 
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(27) 


4.  Base  Pressure  Drag  Increase  Due  to  Presence  of  Fins 


Then*  are  several  primary  iactors  which  determine  the  effect  of  sins 
on  base  pressure.  These  factors  are  fin  location,  thickness  ratio,  aspect  ratio,  profile, 
sweepback  angle,  and  number  of  fins.  Based  on  the  small  amount  of  experimental 
data  available,  it  is  not  possible  to  accurately  account  for  any  of  the  above  factors 
for  a  general  configuration.  However,  order  of  magnitude  effects  of  two  of  the 
variables,  fin  location  and  thickness  ratio,  can  be  estimated  using  References  16 
through  19. 


To  estimate  the  effect  of  fin  thickness  to  chord  ratio  on  base 
pressure,  it  will  be  assumed  the  fins  are  flush  with  the  base.  The  effect  of  the  fins 
not  being  flush  with  the  base  will  be  accounted  for  shortly.  Figure  9  is  a  plot  of 
ACp  ^  A  t/c)  versus  Mach  number.  Here 


AC„ 


“  (CPB\vith  fins  ' 


no  fins 


The  points  above  Mw  =  1.5  were  taken  from  the  data  of  References  10-  and  17 
whereas  those  points  below  =  1.5  were  taken  from  Reference  19.  Again  it 

should  be  emphasized  that  there  are  too  few  data  points  to  put  a  great  deal  of 
confidence  in  this  curve.  For  a  given  Mach  number.  M,,  the  increment  in  base 
pressure  due  to  the  presence  of  fins  is  then: 


(ACn  ),  = 


b  r 


ACn 


L(,/c)JM 


w 


(27) 


The  values  of  t/c  which  (27)  was  derived  for  were  ten  percent  or  less. 

The  work  of  Spahr  and  Dickey17  has  shown  that  if  the  fins 
were  placed  upstieain  a  given  distance  rather  than  flush  with  the  base,  the  effect  on 
base  pressme  is  not  as  great.  Furthermore,  if  the  fins  were  moved  far  enough  from 
the  base,  they  would  have  no  effect  on  base  pressure,  the  amount  of  this  movement 
being  dependent  mainly  on  fin  thickness  to  chord  ratio  and  profile.  As  seen  in 
Figure  10,  this  distance  varies  linearly  with  t/c  up  to  values  of  0.10.  The  curve  in 
the  figure  is  then  extrapolated  from  t/c  =  0.1  to  t/c  =  0.2. 


EXPERIMENT 
A  REF  16 


WITH  BASE 


/ 


FIGURE  10.  DISTANCE  FROM  BASE  WHERE  FINS  DO  NOT  AFFECT 
BASE  PRESSURE. 


Now  if  a  linear  variation  of  (ACp^  )f  is  assumed  between  its  maximum 
when  the  fins  are  flush  with  the  base  and  zero  when  the  fins  are  far  enough  away 
from  the  base,  then  Equation  < 27 }  may  be  modified  in  the  form: 


<ACo  = 


ACp 


B’f  Ki/Ojm.m, 
(ACf.B),.  =  0 


c  e 


t  x 

1  >  0.1  - 


-  <  0.1  - 

c  c 


(28) 


x/e  in  Equation  (28)  is  the  distance  (hi  chord  lengths)  upstream  of  the  base.  This 
empirical  relation  was  derived  only  for  cruciform  fins. 

C  Tail  or  Canard  Alone  Lift 


To  calculate  the  normal  force,  the  wing  is  assumed  to  be  represented  by  a 
flat  plate  of  given  planform  with  zero  thickness  and  camber.  This  assumption  is 
justified  because  thickness  has  only  a  second-order  effect  on  lift  (except  for  thick 
wings  in  transonic  flow)  and  most  missile  configuration0  have  wings  with  zero 
chamber.  The  methods  used  to  compute  the  canard  or  tail  alone  lift  are  different 
for  subsonic,  supersonic,  and  tran-onic  flow  and  will  be  discussed  separate!}  below. 

1.  Subsonic  Flow 

The  basic  equation  of  motion  is  Equation  (1)  with  boundary 
conditions  (3)  and  (4).  The  boundary  condition  (3),  under  the  assumption  of  zero 
thickness  and  chamber.  ma\  be  simplified  to' 


w(x,y)  =  <t>7  -  ~(x 


(29) 


In  addition  to  these  boundary  conditions,  the  Kutta  condition  (which  requires  the 
velocity  on  the  upper  and  lower  surfaces  at  the  trailing  edge  to  be  equal)  is  also 
imposed  for  subsonic  How.  Equation  (1)  may  be  simplified  somewhat  by  using 
(iothert’s  rule' 20)  to  relate  the  compressible  subsonic  normal  force  oi  pitching 
moment  to  the  incompressible  case.  That  is: 


(30) 
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„,< fc./R.a 

(CM  tic. /R. a  = 
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1  -  M2 


(31) 


However,  it  was  assumed  apriori  that  the  only  contribution  to  lift  was  due  to 
angle-of-attaek,  so  Equations  (30)  and  (31)  may  be  simplified  to: 


'CnV„,®,k  ~ 

(C-N  h.PAi  ,pa 

(32) 

i  -  M2 

(33) 

1  -  M2 

The  way  this  rule  is  applied  is  as  follows:  given  a  wing  of  aspect  ratio,  j^R. 
freestrcain  Mach  number,  and  anglc-of-attack,  a ,  the  normal  force  and  pitching 

moment  can  be  obtained  by  calculating  the  normal  force  and  pitching  moment  on 
another  wing  of  aspect  ratio  j3/iR,  Mach  number  zero,  and  anglc-of-attack  ,0a.  Using 
the  above  relations,  the  normal  force  and  pitching  moment  on  a  given  wing  at  any 
subsonic  Mach  number  may  be  found  by  calculating  the  aerodynamics  of  an  affinely 
related  wing  at  zero  Mach  number. 


For  =  0,  Equation  (>)  reduces  to  La  Places  Equation: 

<P  +  <p  +  4>  ~  0 

rx\  ryy 

The  solution  to  Equation  (34)  can  be  shown  to  be:(2l> 


(34) 
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where  ACp  =  Cp^  -  Cp^.  It  is  required  to  determine  the  pressure  loading  ACpover 
the  entire  surface.  Following  Chadwick122',  Equation  (35)  is  first  differented  with 
respect  to  and  the  limit  as  -*•  0  taken.  Tlie  result  is  then  equated  to  the 
boundary  condition.  Equation  (29).  to  obtain: 


a  = 


S7T 


r ACpfXj.y  ,) 

X  -  X, 

]  +  * 

1  (y  -  y,  >2 

ytx  -  x,)2  +  (y-  y,)2_ 

dx,dy, 


(36) 


The  cross  on  the  yj  integral  indicates  a  singularity  at  y  —  y3 ,  in  which  case 
Manglers  principal-value  technique*2 1  *  can  be  applied.  The  details  of  the  solution  of 
the  integral  Equation  (36)  for  AC,.(x.y)  will  not  be  repeated  here  as  they  are  given 
in  detail  in  many  references  (see  for  example.  Reference  22).  Worthy  of  note, 
however,  is  the  fact  that  Equation  (36)  is  an  integral  equation  for  which  the  wing 
loading  ACp  is  to  be  found  as  a  linear  function  of  angle  of  atiack. 

Once  the  span  loading  ACp(x.y)  is  known  over  the  entire  wing 
surface,  the  normal  fore,  at  a  given  spanwise  location  is: 


ACpdx 


The  total  normal  force  for  the  entire  wing  is: 


C 


N 


(37) 


v.’R) 


The  pitching  moment  of  a  given  airfoil  section,  about  the  point  where 
the  wing  leading  edge  intersects  the  body,  is  then  (positive  leading  edge  up): 

XT£ 

xACpdx  (39) 

!•: 
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The  total  pitching  moment  becomes. 


C, 


M 


1 


->ef 


(40) 


If  it  is  desired  to  calculate  the  pitching  moment  about  some  other  reference  point, 
then 


CM  =  C„  +  Cv 
Mo  M  N  g 


*0 


(41 ) 


ref 


where  x0  is  the  distance  from  the  reference  point  to  the  juncture  of  the  wing 
leading  edge  with  the  body.  The  center  of  pressure  of  an  airfoil  section  is: 


x 


ci> 


r 


(42) 


or  of  the  entire  wing 


xci> 


(43) 


Finally,  the  spanwisc  center  of  pressure  of  a  wing  semispan  is: 

b/2  „ 

j  cC^y 

yC5»  =  -jtjj  (44) 

/  tf-.dy 


liquations  (37),  (38),  (39),  (40),  and  (44)  can  be  solved  by  numerical  quadratun.  such  as 
Simpson’s  rule,  with  special  attention  given  to  the  leading  edge  singularity 
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Supersonic  Flow 


Here  again  Fquation  (1)  is  valid  along  with  the  boundary  conditions 
(4)  and  (29).  However,  it  will  be  assumed  in  the  present  work  that  the  wing 

trailing  edge  is  supersonic  (Mach  number  normal  to  wing  trailing  edge  is  greater  than 
one)  so  that  the  Kutta  condition  need  not  be  applied.  Two  cases  will  be 
considered:  supersonic  leading  edges  and  subsonic  leading  cjges.  For  both  cases  the 
How  is  conical  in  nature  from  the  vertex  point  which  is  the  intersection  of  the 
wing  leading  edge  with  the  body. 

a.  Supersonic  Leading  Edge 

Since  the  How  is  conical,  the  flow  properties  are  constant  along 
rays  emanating  from  the  vertex  point  and  lying  behind  the  Mach  line.  Referring  to 
Figure  11  A,  where  the  leading  edge  Mach  line  intersects  the  wing  trailing  edge,  this 
means  that  if  the  flow  properties  at  one  point  on  the  ray  are  known  they  are 
known  all  along  the  ny.  It  is  then  a  matter  of  computing  the  induced  velocities 

ami  hence  iiiting  pressures  on  several  rays  from  point  0.  The  lifting  pressure  on  any 

one  ot  these  rays  in  region  2  is:1 


AC 
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4cv 


"I 


7T 


(45) 


where  17  and  a  were  defined  by  7  A  Ahead  of  the  Mach  line,  in  region  1,  the  flow 
is  two-dimensional  so  the  lifting  pressure  is  constant  at: 


AC 


4a 


(4o) 


If  the  wirg  wore  tapered  to  a  point,  the  above  two  relationships 
would  allow  complete  determination  of  the  lifting  pressures  over  the  wing.  I-'or  most 
practical  cases,  the  wing  tip  *s  not  pointed  so  tip  effects  must  be  accounted  for. 
The  tip  a  feds  the  pressure  within  the  Mach  line  from  the  tip  leading  edge 
(region  3).  Within  this  region,  the  How  is  again  conical  along  rays  from  point  D  so 
the  lilting  pressure  caused  by  the  lip  is  ' 2i) 
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FIGURE  IIA.  FLAT  PLATE  WING  PLANFORM  WITH  SUPERSONIC 
LEADING  AND  TRAILING  EDGES-,  MACH  LINE 
INTERSECTS  WING  TRAILING  EDGE. 
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(47) 


ACp 


3 


, - -tan-ifUsaZ 

TTji  -  V2  V  -0D(1  +  T)) 


Flu1  original  equation  ot  motion.  Equation  ( 1 ),  is  linear  so  that  superposition  of 
Mentions  is  allowed.  Thus,  the  total  lifting  pressure  in  region  3  is: 


AC,,  =  AOP  +  AC,,^ 

It  the  situation  arises  such  that  the  Mach  lines  intersect,  as  shown  at  the  bottom  of 
Figure  11  A.  then  an  additional  pressure  of  AC,.  ?  is  created  in  region  4.  Hence,  the 
total  lifting  pressure  in  region  4  is 


AC,,  -  AC,,  i  +  ACp 


Hie  second  case  to  consider  when  the  leading  edge  is  supersonic 
*'  when  the  leading  edge  Mach  line  intersects  the  tip  as  illustrated  m  Figure  1 1 B. 
the  lilting  pressures  m  regions  1,  2,  3,  and  4  are  calculated  in  the  same  manner  as 
when  the  Mach  line  intersects  the  trailing  edge.  However,  in  region  5  the  lift  to  be 
cancelled  along  the  tip  is  variable  (see  Equation  45)  as  opposed  to  the  constant 
value  cancelled  in  regions  3  and  4  (see  Equation  4b).  This  complicates  the  problem 
'"‘I'  What  in  that  a  summation  (or  integration)  must  now  be  performed  to  calculate 
die  cancellation  lilt  in  region  5.  Referring  to  Figure  111),  the  region  5  is  broken 
down  into  a  finite  number  of  intervals,  in  the  following  manner.  Conical  rays  from 
die  vertex  at  0  are  projected  behind  the  Mach  line  OA  and  intersect  the  wing  tip  at 
equal  intervals  along  the  tip  (for  example,  rays  OJ  and  OK).  Lines  are  then  drawn 
hom  the  points  of  intersection  on  the  wing  tip  parallel  to  the  Mach  line  A1  until 
d'e  point  P(\,v)  is  enclosed.  The  difference  in  ACp  across  one  of  these  intervals  is 
what  must  be  cancelled  throughout  the  region  5.  For  any  given  interval  then,  say 
IK.  this  difference  in  pressuie  is  found  by  applying  Equation  (45)  to  each  of  the 
ia\s  from  0  passing  through  J  and  K.  That  is* 


AC,,  (\,b/2 )  ~  AC„  (\,b/2)  -  ACp  (\.b/2) 
1  h  1 1  1  k 


INTERVAL  JK 


FIGURE  I  IB.  FLAT  PLATE  WING  PLAN  FORM  WITH  SUPERSONIC 
LEADING  AND  TRAILING  EDGES-,  MACH  LINE 
INTERSECTS  WING  TIP. 


ACPjK(x.b/2) 


So _ 

v/1  - 


The  effect  of  this  cancellation  pressure  at  any  point  P(x.y)  is  then: 

AC,,  (x,\ )  =  AC'  (x.b/2)  -  tan  If- V  ------ -  (48) 

5JP  J  K  7T  y-OjpU+ri) 

Now  the  interval  JK  is  any  inten'al  upstream  of  the  Mach  line  from  P  passing 

tluough  the  wing  tip.  If  m  is  the  total  number  of  intervals  upstream  of  P.  then  the 

total  lifting  pressure  coefficient  at  point  P(.\,y)  within  region  5  is: 

I-  U) 

ACP  =  AC,,  +  ACp  +  S  AC,, 

1  2  ‘  4  ...  1  iP 


where  AC.,  is  given  at  each  interval  b\  Hquation  <48). 

The  force  and  moment  coefficients  can  now  be  found  by 
substituting  the  expressions  for  AC,,  in  each  region  into  the  Equations  (37)  through 
(44)  and  performing  the  indicated  operations.  The  integrations  could  be  carried  out 
in  closed  form  for  regions  I  through  4,  but  the  formulas  are  very  lengthy  for  even 
the  simpliest  cases.*’ 2 1  A  more  straightforward  approach  is  to  simply  numerically 
integrate  the  integrals  over  the  entire  surface,  particularly  since  the  pressures  in 
region  5  lure  to  be  numerically  integrated. 

b.  Subsonic  Leading  Edge 

For  subsonic  leading  edges,  the  velocity  and  lifting  pressures 
approach  infinite  values  near  the  leading  edge.  The  solution  for  ACp  is  complicated 
somewhat  by  this  singularity,  but  due  to  the  conical  nature  of  the  flow  from  the 
urUx  points  0  and  t)  (Figure  12).  a  closed  form  solution  can  he  obtained.  For 
region  I  of  Figure  12  this  is:*2 1  ) 


4a- 


(49) 


AC,, 


E  tan  AiN/T1^ 


where  E  is  a  constant  dependent  only  on  the  leading  edge  sweep  angle  and  Mach 
number. 

That  is. 


E(k»  = 


f J 

-'n 


J 1  -  k2  sin2  OdO 


(50) 


and  k  =  J\  -  0~Uc.\\~  Al(.  .  The  value  of  the  complete  elliptic  integral  H(k)  has 
been  tabulated  and  appears  in  standard  mathematical  handbooks. 


The  lift  to  be  cancelled  at  the  tip  is  variable  as  was  the  case  for 
supersonic  leading  edges  when  the  Mach  line  intersected  the  tip.  However,  in  this 
case  the  integration  can  be  carried  out  in  closed  form.  Thus,  the  cancellation  lifting 
pressure  is:'-3’ 


v/here 
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{ 2j8(b/2-y) 

W/2e  ft'1) 

K(k^ 
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(51) 
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Here  F(0,.(?.,)  and  E(0,,0.,)  are  elliptic  integrals  of  the  first  and  second  kind, 
respectively,  which  again  can  be  found  in  any  standard  mathematical  handbook.  The 
complete  elliptic  integral  K  is  related  to  F (0,,02)  by  K  =  F(0,,7r/2)  and  the 
complete  elliptic  integial  E  of  Equation  (50)  is  related  to  E(0,,0.,)  by 
E  =  E(0,  ,tt/2). 

The  total  lifting  pressure  coefficient  at  any  potnt  within  region  2 

is  then 


AC,,  =  AC,,  ^  +ACp2 

The  integrations  for  spanwise  lift  and  pitching  moment  can  be 
obtained  by  integrating  Equations  (37)  and  (39)  numerically.  However,  caution  must 
be  taken  in  these  integrations  because  for  subsonic  leading  edges,  the  lifting  pressure 
goes  like  l/yHT  near  the  leading  edge.  In  this  vicinity,  more  mesh  points  must  be 
added  to  the  flow  field  to  get  an  accurate  integration. 

The  total  wing  normal  force  and  pitching  moment  can  be  easily 
obtained  by  numerical  integration  of  Equations  (38)  and  (40).  The  cho.'dwise  center 
of  pressure  is  then  found  from  Equation  (43)  and  the  spanwise  center  of  pressure 
from  Equation  (44). 

3.  Transonic  Flow 

As  mentioned  earlier,  airfoil  thickness  has  a  second  order  ‘fleet  on 
lift  in  subsonic  and  supersonic  flow.  However,  this  is  not  true  in  transonic  flow  so 
the  assumption  of  a  flat  plate  with  zero  thickness  is  no  longer  valid  except  'or  very 
thin  wings.  Furthermore,  as  approaches  unity.  Equation  (1)  cannot  be  applied 
because  the  term  ( 1  -  M2  )<t>  ,  becomes  of  the  same  order  as  nonlinear  terms  which 
have  been  neglected  in  deriving  this  equation  from  small  perturbation  theory. 

Recent  progress  iri  the  field  of  transonic  aerodynamics  has  greatly 
advanced  the  state-of-the-art.  However,  at  present,  practical  methods  for  flow  field 
computation  are  still  severely  limited.  For  example,  solutions  for  three-dimensional 
swept  and  tapered  wings  with  thickness  do  not  exist,  even  in  approximate  form.  In 
light  of  these  considerations,  an  empirical  approach  to  wing  lift  in  transonic  flow 
will  be  used.  The  method  adopted  is  that  presented  in  DATCOM.124'  This  procedure 
accounts,  in  an  empirical  manner,  for  sweep.  Mach  number,  aspect  ratio,  and 
thickness  ratio  but  not  for  airfoil  section. 
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To  apply  the  above  empirical  procedure,  the  force  break  Mach  number 
is  found  from  Figure  i3A  tor  a  wing  of  zero  sweep  and  corrected  for  sweep  by 
Figure  13B.  The  lift  curve  slope  at  the  force  break  Mach  number  is  then  computed 
by  a  simple  expression  derived  from  lifting  line  theory:* 251 


(C 


a  fb 


_ 2jt/R _ 

2+  [/R2(02  +  tan2  AI/l)  +  4]’/2 


(59) 


This  value  is  corrected  to  agree  with  experiment  according  to  Figure  14A.  The 
abrupt  decrease  in  lift  curve  slope  associated  with  thick  wings  is  approximated  by 
the  relation: 


where  a/c  is  given  in  Figure  14B.  The  Mach  number  at  point  a  is: 

Ma  =  Mfb  +  0.07 

The  subsequent  rise  in  CNq  to  a  value  at  point  b  is 

<cN  1  -  n-b/o(cN 

(XV 

where  b/c  is  given  in  Figure  140.  The  Mach  number  at  point  b  is 


(60A ) 


(61 ) 


Mb  =  Mfb  +0.14 


(61  A) 


The  normal  force  curve  slope  at  M  >  1.2  is  calculated  by  supersonic  it. in  wing 
theory  and  for  M  <  Mfb  by  lifting  surface  theory.  From  this  empirical  correlation 
one  obtains  (CNft)  ,  (Cj^K  (CNft)()  at  the  Mach  numbers  Mfb,  Mfb+.07,  and 
M,  +  .14.  For  values  of  Cm  in  transonic  (low  in  between  these  Mach  numbers. 

In  <x 

interpolation  is  used. 


THICKNESS  RATIO (% CHORD) 


FIGURE  I3A.  FORCE- BREAK  MACH  NUMBER  FOR  ZERO  SWEEP 


FIGURE  138.  SWEEP  CORRECTION  FOR  FORCE-BREAK  MACH  NUMBER 


FIGURE  13.  TRANSONIC  FORCE-BREAK  MACH  NUMBER 
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THICKNESS  RATIO  (%  CHORD) 

FIGURE  I4A.  CORRECTION  TO  LIFT-CURVE  SLOPE  AT  FORCE-BREAK  MACH 
NUMBER 
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THICKNESS  RATIO  (%  CHORD) 

FIGURE  I4B.  CHART  FOR  DETERMINING  LIFT- CURVE  SLOPE  AT 
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THICKNESS  RATIO  (%  CHORD) 

FIGURE  14 C.  CHART  FOR  DETERMINING  LIFT-CURVE  SLOPE  AT 

FIGURE  14.  CHART  FOR  DETERMINING  TRANSONIC  LIFT-CURVE 
SLOPE  AT  Ma  AND  Mb 


The  center  of  pressure  of  wing  alone  lift  is  usually  around  the  quarter 
chord  for  subsonic  How  and  half  chord  for  supersonic  flow.  Transition  front 
subsonic  to  supersonic  flow  is  assumed  to  occur  in  a  linear  fashion  between  values 
calculated  at  slightly  less  than  M)b  and  >  1.2.  The  pitching  moment 
coefficient  derivative  of  the  wing  alone  is  then: 


'M 


~  xo  p  Cn  a 


(62) 


D.  Interference  Lift 

Interference  lift  is  broken  down  into  three  parts:  lift  of  canard  or  tail 
due  to  presence  of  body,  lift  of  body  due  to  presence  of  wing  or  canard,  and 
vortex  lift  on  tail  due  to  wing  shed  or  body  shed  vortices.  The  methods  used  to 
calculate  the  interference  lift  components  are  essentially  those  presented  in 
Reference  26.  The  necessary  equations  for  the  various  calculations  will  be  repeated 
herein,  but  for  the  details  of  the  derivations,  the  reader  is  referred  to  the  above 
reference. 


The  method  used  by  iVlorikawa(27)  for  presenting  lift  interference  is 
adopted.  He  defines  the  wing  alone  as  the  exposed  half  wings  joined  together.  The 
lift  of  the  combination  is  related  to  the  lift  of  the  wing  alone  by  the  factor  K(.  so 
that 


Lc  Kc  Lw 


(63) 


Kc  is  actually  composed  of  three  components. 
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Kw(.I)  +  KB 


(64) 


which  are  the  ratios  of  the  body  lift  in  presence  of  the  wing,  wing  lift  in  presence 
of  body,  and  body  lift  to  the  lift  of  the  wing  alone.  That  is: 
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We  will  not  be  concerned  with  KB  as  it  was  found  in  Reference  i  as  discussed 
earlier.  The  factors  Ku(w)  and  Kw(B)  are  found  from  a  straightforward  application 
of  slender  body  tneory  and  are: 

(1  -  r2/s  2 )  -  2/tr  -|(1  +r4/s4)[V2  tan~ 1  'A(s/r-  r/s)  +  tt/4| 

Kb(w)  =  (l  -  r/s 


r2/s2[(s/r-r/s)  +  2  tan"  1  (r/s))} 
(1  -  r/s)2 


(66) 
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(67) 


r2/s2((s/r-  r/s)  +  2  tan"1  (r/s )!} 
( 1  -  r/s  )2 


Kw([j)  anc^  ^E(w)  Equations  (66)  and  (67)  are  related  to  each  other  by: 


^w(It)  +  ^B(w)  “  ^  +r/s 


(68) 
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Values  of  these  parameters  are  plotted  in  Figure  15.  Equation  (66)  is  used  to 
calculate  K„,  ,  for  all  subsonic  Mach  numbers  and  for  supersonic  Mach  numbers  il 
Ml(l  +  X)  U/rn/3)  +  1 1  <  4.  If  /3/R(l  +X)Il/(mp)+  l]  >  4,  Nielsen'-8 >  has  shown 
that  for  no  afterbody  (sec  Figure  16),  a  more  appropriate  value  of  K„(w)  can  be 
found  from: 


8  /ftf 

kb(.,i«c»Pw"x  +i)ls,r’ " =  U 


cos" 1 


m/3  +  cr  /(/3d) 


+  mcr  Id 
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P  \/3d  J  \mp) 


+  m/3  FM  \J m2/32  -  1  shf 1  f— j  -  cos  If  1  ^ } 


(69  A) 


when  j3m  >  1  and  cr//3  >  d.  If  pm  <  1  and  er/p  >  d,  K„(w)  is  found  from: 


lbv^m/T-  //3d\  1/  ,  111cr  \  //£ r_  +1 

"Is'r- " =  ^TTit)  v +  ”i  Mud  )v  d 


(69  B) 


A  (m/5)3/2  +  v/3m  +  1 ) 
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\T^ 


If  c  //3  <  d  in  Equations  (69A)  or  (69B),  cf//3  is  set  equal  to  d.  For  an  infinite 
afterbody  behind  the  wing  or  canard  (see  Figure  16),  Equations  (69 A)  and  (6913) 
are  replaced  by: 
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SHORT  AFTERBODY 


INFINITE  AFTERBODY  ( XAFT  =  d  COT  M) 


NO  AFTERBODY  (  XApT*  0) 


FIGURE  18,  DETERMINATION  OF  Kg(w)  FOR  HIGH -ASPECT-RATIO 
RANGE  AT  SUPERSONIC  SPEEDS. 
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for  ni0  >  1  and  by 
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for  in0  <  1  For  cases  in  between  no  afterbody  and  an  infinite  aferbody.  KR(wJ  maybe 
found  approximately  by  linearly  interpolating  between  the  values  given  by  Equations  (ft1) ) 
and  (70).  That  is,  referring  to  Figure  16, 
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Mere  the  subscripts  SA,  IA,  and  NA  refer  to  short,  infinite,  and  no  afterbody, 
respectively. 

In  addition  to  the  normal  force  due  to  angle  ol  attack,  there  is  a 
normal  force  created  due  to  wing  deflection  5.  The  lift  of  a  wing  due  to  a 
deflection  6  is  given  by: 


cs  =  kWIB1(CK  >  5 

NwlB>  No  % 


1  1 1 1  <Wr  +  1  >:  +  +  H 2  ,  f  (s/r)2  -  ll 

»■  (4  (s/ri-  (s/r)2(Wr-UJ  [_(x  /r>2  +  IJ 

-.1^11?+  ingliilL  L-<  ri^i^  J 

vr(>/r  *l)  (>/r)i(s  /r  -  1  )2  y  f(»/r)J+lJ  J 

.  iL'iLUisin-i  p^nL_Ll  +  _JL_l0ffe£)LLi]j 

</rls  /r+l »  Lf,»/r>2  U  fa/  r-  I)2  L  2s/r  JJ 


!  ik1  hit  on  l he  hodv  due  to  the  control  deflection  is 


NBlw,Tki<»^N®>w  5 


w  lie  re 


kW(ln  ;unl  kI)fw)  are  given  in  Figure  15  as  a  function  of  r/s. 

Stnctly  speaking.  Equations  (66).  (67).  and  (73)  apply  to  slender 
winp-l'oih  .on figurations  which  have  low  aspect  ratio  wings.  Furthermore,  the 
trailing  edges  of  these  wings  must  be  perpendicualr  to  the  body  axis  or  swept 
forward  hut  cam  ot  be  swept  back.  Ncilscn(26)  has  shown  that  good  results  can  be 
obtained  tor  the  interference  lift  components  for  high  aspect  ratio  wings  even 
though  the  slender  body  approximation  is  made.  However,  if  the  wing  trailing  edge 
is  swept  back  the  previous  slender  body  formulas  arc  not  applicable  without  further 
stmlv.  For  wings  which  have  swept  back  trailing  edges,  and  where  slender  body 
tneorv  jv  used  to  calculate  interference  lift.,  it  seems  reasonable  to  assume  the 
incmavd  lilt  op  the  wing  due  to  the  presence  of  the  body  is  distributed  evenly 
along  the  wing  chord  and  is  concentrated  near  the  root.  With  this  assumption,  the 
i :  *  -if.  lift  components  of  a  swept  back  wing  may  be  approximated  by: 
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Best  Available  Copy 


K,bj]u  =  H  (K,b,],-,-)G 

[kw(B)]jI  ~  1  +  (j^wUJ)]  I  _1)G 

[kI!  <  w  »]  1 1  "'([kw(inJ[  -  [k  w  ( B  )]  l)  ^ 


where  G  -  (<-'r),,/(cf),  and  the  subscript  Ii  refers  to  the  actual  wing  being 
considered  and  subscript  I  refers  to  the  wing  for  which  the  slender  body  liquations 
(60),  (67),  and  (73)  assumes.  Figure  17  indicates  how  the  above  procedure  is  carried 
out. 


The  negative  lift  caused  by  dow awash  of  the  canard  shed  vortex  on  the 
tail  is  given  by: 


IKw(B)  sin  a  +kw(lj)  sin  S|,(s  -  r)wSw 
2»*w(fw  -  rw  >Srof 


(7oi 


This  equation  is  determined  from  line-vortex  theory  assuming  only  one  trailing 
vortex  per  forward  wing  panel.  The  lateral  location  of  the  vortex,  f  ,  measured 
from  the  body  center  line  is: 


‘w  ~  '  +  2(  (C*  c  VI, 


(77) 


I  he  span  loading  at  the  root  (C'ncir  and  normal  force  coefficient  are  known  for  the 
particular  wing  geometry  as  discussed  previously.  The  interference  factor,  i.  of 
1  quation  (70)  is. 


4b 


DESIRE: 


FIGURE  '7A.  WING  FOR  WHICH  INTERFERENCE  LIFT  IS  DESIRED 


FIGURE  178.  ASSUMED  SLENDER  BODY  REPRESENTATION 


FIGURE  17.  PROCEDURE  USED  TO  CALCULATE  INTERFERENCE 
LIFT  FOR  WINGS  WITH  SWEPT  BACK  TRAILING 
EDGES  WHEN  SLENDER  BODY  THEORY  IS  USED. 
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where 


L[X, = 
l  b  b  b 


s-  rX  -  '( !  -  A)  /lr+(f-s)2 
-fin 


2(s  -  r) 


Ur+(f-  r)2 


I  -  X 


( s  -  r)  +  h  tan”  ! 


f-s 


ii  tan 


-i  A- ^ 
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■  IL-  J!l_ 
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and 


h  =  -  3(er)xv/4sm  6  +  |CT  -  <fw  -  3(cr)w/4|sino; 

Here  it  is  assumed  that  the  hinge  line  of  tiic  wing  is  at  the  quarter  chord. 

The  only  remaining  lift  is  the  negative  lift  on  the  tail  due  to  the  body 
shed  vortex.  This  is 


-N|l(v> 


-4T1 

«wV- 


r2 


_  r  +  —ll - 

T 


(79) 
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where 


r_ 

V 


=  k*(P)5  (r 

4(f  -r  )  ‘  N< 

'  W  U.*  ' 


K  Stt. 


The  center  of  pressure  of  the  wing  lift  in  presence  of  the  body,  canard 
shed  vortex  lift,  and  body  shed  vortex  lift  are  all  assumed  to  act  at  the  center  of 
pressure  of  the  wing  alone.  The  center  of  pressure  of  the  body  lift  caused  by  the 
wing  is  Mach  number  deper  nt.  Complicated  expressions  could  be  derived  for 
center  of  pressure  of  this  rift  component  in  supersonic  speeds  analogous  to 
Equations  (69)  and  (70).  However,  a  good  engineering  approximation  is  to  assume 
that  the  center  of  pressure  is  at  the  centroid  of  the  cross-hatched  area  in  Figure  16. 
This  eliminates  the  need  for  considering  separately  the  afterbody,  no  afterbody,  and 
short  afterbody  cases.  Thus,  for  supersonic  flow,  the  center  of  pressure  of  the  body 
lift  caused  by  the  wings  is  (measured  from  the  wing  leading  edge): 


c?d/2-j32d3/6+  xAFT(2d-  xAFT/P)x-,/ 2 
dcr  -/3d2/2  +  xAFT(2d-  xAI.T/0)/2 


where 


d/2-  xAF-F/(3/3) 
-  d  '  x A  FT 


(80A) 


for 


/3d  -  xAFT  >  0  and  t/3d  -  cf  -  xAFT)  <  0. 

If  /3d-  xAFT  >  0  and  (/3d-  cf-  xAFT)  >  0,  then: 

x  =  Il2JS>  likElf  -  <cr  +  2/3xapt)x2aft 
Cp  ^Cr  +  XA  FT  ^  "  XA  FT 


(80B) 
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Finally.  if  (5d  *  xA|1  <  0.  then 


c  r  +  ffd 


CP 


tson 


Fquatioa  (80A)  represents  the  ease  where  the  Mach  line  from  the  trailing  edge 
intersect  the  base  v.ncrvas  that  from  the  leading  edge  intersects  the  opposite  side. 
Iquation  (SOB)  is  used  when  both  leading  and  trailing  edge  Mach  lines  from  the 
wing  root  chord  intersect  the  base  whereas  tquation  (800  is  used  when  both  Mach 
lines  intersect  the  opposite  side. 

For  subsonic  and  transonic  How,  there  is  much  less  tendency  for  the  lift 
to  be  carried  aft.  For  subsonic  flow  it  is  reasonable  to  assume  the  body  carry-over 
lift  acts  at  the  center  of  pressure  of  the  wing  alone  lift.  In  transonic  flow,  the 
center  of  pressuie  is  assumed  to  vary  linearly  between  its  value  at  the  force  break 
Mach  number  and  that  computed  for  M  =  1.2  by  Fquation  (80). 

Fi.  Summary  Configuration 

The  total  normal  force  of  the  entire  configuration  is: 

^  N  ~  +  [OXv(B)  +  a+  (^wrm  +  ^B(w))  f(-Na\v 

(81 ) 

1  [(KT(B)  +  Kii!T))a](CNu)r  +cn1(V)  +  cnIMV) 

Petining 


C‘\v,B)  '  [KwtB)a  +  kw(B)5](CNjw 
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B(w) 


fiw  +  ^cp^B(w) 


*T(B)  ”  ^T+^Xcp^T 
^B(T)  =  +^Xcp^B(T) 


the  center  of  pressure  referenced  to  the  nose  tip  may  be  written  as: 


C  "  |  eBCNB+CNw(B)ew(B)+CNB(w)CB(w)+CNT(B)CT(B) 
+  CNb(t)CB(T)  +  (CNT(V)+CNB(V))1!t1/CN 


(82) 


The  pitching  moment  about  the  nose  is  then: 


(83) 


A  summary  of  the  various  theoretical  and  empirical  procedures  used  to 
calculate  the  static  aero '  "f  the  wing  or  tail  alone  and  interference  effects 

is  given  in  Figure  t  .oady  indicated,  a  summary  of  the  body  alone 

procedures  was  given  .  .0are  2. 
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Ill.  RESULTS  AND  DISCUSSION 


A.  Comparison  With  Exact  Linear  Theory 
1.  Wing  Wave  Drag 

For  wings  which  have  simple  planforms,  exact  linear  theory  solutions 
exist  by  whidi  the  present  numerical  solutions  may  be  compared.  Simple  planforms 
are  defined  as  wings  with  double  wedge  or  biconvex  airfoil  sections  with  the  same 
airfoil  svaion  and  thickness  ratio  all  along  the  span.  Simple  planforms  also  have 
sharp  leading  and  trailing  edges. 

Figure  19  compares  the  pressure  coefficient  of  the  present  method 
with  that  of  exact  linear  theory'  for  a  biconvex  airfoil  design/131  The  pressure 
coefficient  has  been  divided  by  t/c,  the  thickness  to  chord  ratio,  since  it  is  constant 
all  along  the  span.  As  seen  from  the  figure,  the  numerical  solution  gives  essentially 

exact  results  when  compared  with  the  solution  of  Reference  i3. 

Section  wave  drag  along  the  span  of  a  wing  with  a  double  wedge 
airfoil  section  is  given  in  Figures  20  and  21.  Figure  20  is  for  a  wing  tapered  to  a 

point  (no  tip  effects)  and  Figure  21  is  for  a  wing  with  fifty  percent  taper.  The 

exact  analytical  solution  is  taken  from  Reference  31.  Section  wave  drag  at  seventeen 
spanwisc  stations  were  computed  but,  as  will  be  discussed  later,  only  nine  are 
necessary  for  reasonably  accurate  total  wing  drag.  Again  confidence  in  the  numerical 
solution  is  gained  by  the  near  perfect  agreement  with  analytical  solutions. 

The  final  two  figures  comparing  numerical  drag  calculations  of  simple 
planforms  with  :\act  solutions  give  total  wing  wave  drag  for  biconvex  and  double 
wedge  airfoil  sections.  'Flic  analytical  solutions  come  directly  from  the  charts  of 
Reference  30  Figure  22  is  for  the  double  wedge  airfoil  wing  and  Figure  23  for  the 
biconvex  airfoil  section  wing.  Both  wings  have  zero  taper.  Slight  discrepancies  of  up 
to  three  to  five  percent  between  the  numerical  and  analytical  solutions  exist  at 
some  Mach  numbers  as  seen  from  the  figures.  This  is  due  to  the  numerical 
integration  which  smoothes  out  kinks  in  the  pressure  curves  and  due  to  truncation 
errors.  However,  in  light  of  the  accuracy  of  the  exact  linear  theory  compared  with 
experiment,  this  discrepancy  is  quite  acceptable  for  design. 

All  of  the  preceding  cases  were  for  sharp  leading  and  trailing  edge 
wings  witli  simple  planforms.  For  more  complicated  configurations,  closed  form 
linear  theory  solutions  for  drag  have  not  been  obtained.  However,  the  present 
method,  which  determines  the  source  strength  locally  can  still  be  applied.  For 
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FIGURE  19.  PRESSURE  DISTRIBUTION  ALONG  THE  SPAN  OF 
A  SWEPTBACK  WING  WITH  A  SYMMETRICAL 
CIRCULAR  ARC  AIRFOIL;  XsSO®,  = 

Ct/  Cf  SL 
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FIGURE  22.  WAV 
AIR 


example,  consider  the  wing  shown  in  Figure  24.  Tins  wing  has  3  variable  airfoil 
section  and  variable  thicknc"  all  along  the  span.  Its  design  is  dictated  by  structural, 
rather  than  aerodynamic,  considerations.  Also  shown  on  the  figure  is  a  plot  of 
pressure  coefficient  along  the  chord  lor  two  spanwise  locations. 

A  second  example  of  wings  not  covered  in  exact  linear  theory  is  the 
wing  with  blunt  leading  and  trailing  edges  Figure  25  gives  the  total  drag  coefficient 
(wave  plus  trie  turn  plus  trailing  edge  separation)  for  a  biconvex  airfoil  with  blunt 
leading  and  trailing  edges  and  for  three  Mach  numbers.  Although  no  experimental 
data  exist  to  compare  the  theory  with,  the  trends  are  generally  what  one  would 
*peet.  That  is  bluntnexx  has  a  large  elfcct  on  drag  at  high  supersonic  Mach 
numbers  but  at  low  supersonic  Mach  numbers,  the  effect  is  not  so  large.  However, 
at  these  low  supersonic  Mach  numbers  the  linear  theory  results  for  very  small 

bluntness  ratios  are  questionable.  This  is  because  in  the  neighborhood  of  the  leading 
edg.  where  the  slope  is  near  the  maximum  of  fifteen  degrees  (linear  theory  starts  at 
slopes  of  fifteen  degrees  where  Newtonian  theory  ends),  the  pressure  coefficients 
predicted  In  linear  theory  arc  quite  large.  The  discontinuity  between  the  pressme 
eoelficient  ot  Newtonian  theory  and  that  of  perturbation  theory  is  thus  much  larger. 
However,  as  the  blindness  ratio  is  increased  the  area  over  which  linear  theory  is 
applied  decreases  and  hence  one  expects  the  results  to  be  somewhat  better. 

2.  Wing  Lift 

Three  eas-.w  are  considered  as  test  cases  to  compare  the  numerical 

solutions  with  closed  form  analytical  solutions  such  as  presented  in  References  22 
and  23.  These  include  a  wing  with  subsonic  leading  and  supersonic  trailing  edges,  a 
wing  with  supersonic  leading  and  trailing  edges  with  the  Mach  line  intersecting  the 
tip.  and  a  wing  with  supersonic  leading  and  trailing  edges  with  the  Mach  line 
intersecting  the  trailing  edge.  Hach  of  these  cases  is  sufficiently  different  so  as  to 
cheek  all  methodology  developed  in  the  numerical  calculation  of  lift  on  wings  in 
supersonic  flow.  The  theory  used  in  subsonic  flow  has  been  previously  verified  by 
Chadwick. '■  -  ’  Test  data  must  suffice  to  determine  the  correctness  and  accuracy  of 
the  transonic  methodology. 

Figure  2b  compares  the  wing  loading  along  the  span  ol  a  wing  with 
subsonic  leading  an. I  supersonic  trailing  edges  corresponding  to  the  first  lest  case 
above.  The  Mach  lines  and  regions  where  linear  superposition  of  solutions  occur  arc 
indicated  on  the  fimire.  As  is  apparent,  the  numerical  solution  is  the  same  as  the 

analytical  solution  to  the  accuracy  .»1  nlotting  the  data. 

hguies  27  and  2b  compare  the  analytical  and  numerical  span  loading 
calculations  lor  a  so  per sonic  leading  and  trailing  edge  wing  In  Figure  27.  the  Mach 
line  just  barely  intersect-*  the  Hailing  edge  whereas  in  Figure  28  the  Mach  line 
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NUMERICAL 


riGURE  26.  WINS  WITH  SUBSONIC  LEADING  AND  SUPERSONIC  TRAILING  EDGE 


ANALYTICAL 

NUMERICAL 


FIGURE  28.  WING  WITH  SUPERSONIC  LEADING  AND  TRAILING  EDGES  ;  MACH  LINE 
INTERCEPTS  TIP 


intersects  the  tip.  !n  botli  cases,  the  analytical  solution  is  reproduced  by  the 
numerical  procedure. 

B.  Comparison  With  Experiment 

Several  example  casts  are  chosen  to  compare  evith  experimental  data 
so  as  to  place  some  degree  of  accuracy  on  the  theory.  The  first  of  these  is  a  ten 
caliber  missile  with  clipped  delta  tail  fins.  'Hie  experimental  data  are  taken  from 
Reference  34,  which  gives  the  static  aerodynamics  for  0.8  <  Mm  <  1.3.  Figure  29 
compares  theoretical  drag  coefficient,  normal  force  coefficient  derivative,  and  center 
of  pressure  with  the  experiment  as  a  function  of  Mach  number  and  for  a  =  1°. 
Recall  that  for  M„  >  1.2.  the  lift  and  drag  (except  for  base  drag)  was  calculated 
numerically  whereas  for  0.8  <  M„  <  1 .2,  the  theory  consists  of  mostly  empirical 
procedures.  For  <  0.8,  the  wing  lift  is  calculated  analytically  but  most  other 
force  components  arc  computed  empirically.  With  the  exception  of  the  normal  force 
coefficient  slope  at  =  0.8  and  0.85.  the  theory  is  well  within  ten  percent  of 

experiment.  The  maximum  error  in  center  of  pressure  for  this  configuratbn  is  five 

percent  of  the  length  or  half  a  caliber. 

The  next  case  considered  is  again  raken  from  Reference  34  where  the 
same  body  geometry  was  tested  with  several  different  fins  present.  The  case 
considered  is  for  aspect  ratio  two  rectangular  fins  flush  with  the  base.  The  body  is 
the  same  as  that  shown  in  Figure  29.  Figure  30A  presents  the  small  angle-of-attack 
results  for  Cn.CNft,  and  xc()  as  a  function  of  Mach  number  and  Figure  30B  gives 
C|) ,  Cv.  and  x  as  a  function  of  oc  for  Mach  number  1.3.  In  both  Figures  30A 
and  30B.  the  results  are  acceptable  from  an  accuracy  standpoint,  although  the  diag 
coefficients  at  M^  =  1.20  and  1.3  are  off  more  than  for  most  configurations 
considered.  The  reason  most  likely  lies  with  the  empirical  estimation  of  base 
pressure  increase  due  to  presence  of  fins.  At  M,„  =  1.3,  this  component  accounts 

for  about  0.05  of  the  drag  coefficient,  which,  for  this  wing-body  configuration, 

appears  high. 


The  final  example  chosen  is  a  complex  canard-body-tai!  configuration. 
The  body  nose  is  sixty  percent  blunt  with  two  ogive  segments  and  a  0.7  caliber 
boattail.  The  canard  has  an  aspect  ratio  of  two  with  a  swcepback  angle  of  1 5°.  Its 
shape  consists  of  a  sharp  wedge  leading  edge  with  a  constant  thickness  section  following. 
The  trailing  edge  is  truncated  parallel  to  the  leading  edge.  The  tail  has  an  aspect 
ratio  of  four  with  cylindrical  leading  and  trailing  edges  and  where  A,  =  30°, 
A,  =  22.5°,  A3  =  37°,  and  A4  =  30°.  The  tail  thickness  to  chord  ratio  also  vanes 
along  the  span.  The  detailed  canard  and  wing  geometry  listed  above  is  not  needed 
in  calculating  lift,  but  it  must  be  known  for  drag  computations.  The  results  of  the 
calculations  for  this  configuration  are  shown  in  Figure  31.  Figure  3’ a  gi\es  the 
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FIGURE  29.  STATIC  AERODYNAMICS  OF  A  MISSILE  CONFIGURATION 
1.O, J_,  =53.1°,  X  s0.5  ,oc  =  1° 
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?ODYNAMICS  OF  A  MISSILE 
HON;  ^=2.0,X  *0°,  A*  1.0 


ANGLE  OF  ATTACK  -  CC  (DEGREES) 

FIGURE  30B.  STATIC  AERODYNAMICS  OF  A  MISSILE 
CONFIGURATION-,  /R=2.0t X=0°,  A  =10, 


normal  force  and  center  of  pressure  for  M„.  =  1.6  and  at  various  angles  of  attack. 
Four  curves  are  shown  in  the  figure.  canard-body -tail  with  canards  deflected  up  In 
ten  degrees,  canard-bodv-tail  with  no  canard  deflection,  body-tail,  and  finally  body 
alone.  Several  points  are  worthy  of  note  in  this  figure.  First  of  all  the  body  alone 
solution  agrees  very  well  with  the  unpublished  experimental  data  up  to  a  =  In 
Above  or  =  1  o',  the  theory  is  low  which  is  probably  due  to  not  tahine  into 

account  Reynold's  number  effect  in  the  body  crossflow  drag  coefficient.  The  next 
point  is  that  for  this  configuration,  the  tail  lift  is  about  ten  percent  loo  high  and 

the  canard  lift  about  15'?  too  low  so  that  the  total  lift  agrees  almost  perfectly 

with  the  expeiimental  data  up  to  the  point  where  stall  begins  to  occur  (o  A-’  14°). 

This  in  turn  causes  the  center  of  pressure  to  be  more  rearward  than  the 

experimental  data  suggest  by  about  half  a  caliber.  It  is  suspected  that  the  theory 
being  high  for  the  high  aspect  ratio  tail  and  low  for  the  moderate  aspect  ratio 

canard  is  due  to  the  flow  field  interaction  effects  from  the  complex  configuration 

and  will  not  in  general  be  true  for  other  cases.  However,  it  docs  indicate  that  the 

theory  can  he  used  quite  effectively  in  design,  even  for  quite  complex  wing-bodv-tail 
geometries.  The  final  point  to  be  emphasized  from  Figure  3 1 A  is  the  fact  that  no 
attempt  has  been  made  to  predict  stall  characteristics.  As  seen  in  the  figure,  for  this 
configuration,  stall  occurs  around  a=  15°  at  -  1.6.  However,  if  the  wing 
thickness  or  ficestream  Mach  number  is  changed  the  stalling  angle  of  attack  will  also 
change. 


The  drag  characteristics  for  this  same  missile  are  shown  m 

Figure  3 IB.  The  drag  is  shown  as  a  function  of  Mach  number  and  again  the  total 
force  is  broken  down  into  its  components;  body  alone,  body  tail,  and 

canard-body-tail  The  body  alone  drag  is  acceptable  in  supersonic  and  subsonic  flow 
but  is  unacceptable  in  transonic  flow  where  the  empirical  nature  of  the  theory  do<.s 
not  account  for  nose  bluntness  correctly.  The  wing  alone  drag  shown  at  the  bottom 
figure,  includes  the  increase  in  base  drag  clue  to  tails.  This  causes  the  tail  drag  to 
be  high  because  the  theory  predicts  this  base  drag  increase  to  be  significantly  higher 
than  the  experimental  data  suggest.  However,  the  body-tail  drag  is  still  within  flu* 
ilOf?  catagory.  Finally,  the  canard  drag  shown  at  the  bottom  figme.  is  added  to 
the  body-tail  drag  and  the  overprediction  of  tail  drag  is  compensated  omewh.it  by 
the  under  prediction  of  canard  drag. 

As  was  mentioned  earlier,  the  tail  has  a  cylindrical  leading  edge  lot 

which  the  combined  Newtonian  perturbation  theory  must  be  used  to  calculate  the 

pressure  coefficient  Figure  310  presents  this  pressure  coefficient  in  the  vicinity  <>| 
the  leading  edge  lor  Mach  number  two  and  at  a  spanwise  station  ol 
y/(b/2)  =  12.5%.  The  discontinuity  in  pressure  coefficient  caused  In  the  diflereiue 
in  Newtonian  and  perturbation  theory  estimates  occurs  at  a  nu'.ch  point  <>t 
\/c  ■  0.000. 
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FIGURE  3IA.  NORMAL  FORCES  AND  CENTER  OF  PRESSURE  OF 
A  MISSILE  CONFIGURATION;  ,(R  =2  M  =1  S 

T  ®  '  OO  * 


FIGURE  31 B.  DRAG  OF  A  MISSILE  CONFIGURATION  AND 
ITS  COMPONENTS 


y/(b/2)  =  12.5% 


FIGURE  31C.  PRESSURE  COEFFICIENT  ON  LEADING  EDGE  OF 
A  BLUNT  LEADING  EDGE  WING  j  JL|  =30°,  Moo  =  2, 
b a 0.92', CR=  .23',  R^.0021',  tR  =.0167*  tf  =  .0042' 
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C.  Computational  Time  and  Cost 

Although  the  method  presented  herein  for  computing  aerodynamics  of 
guided  or  unguided  projectiles  actually  consists  of  several  rather  complicated 
theoretical  and  empirical  procedures,  the  cost  to  obtain  force  and  moment 
predictions  is  relatively  small.  For  example,  the  most  complicated  configuration 
considered  to  date  was  the  canard-body-tail  shape  in  Figure  31.  To  compute  the 
aerodynamics  at  ten  Mach  numbers  for  a  smali  angle-of-attack  takes  about  five 
minutes  on  the  CDC  6700  compute*  or  costs  about  $75.  Aerodynamics  at 
supersonic  Mach  numbers  costs  about  twice  as  much  per  Mach  number  than  ;•! 
subsonic  or  transonic  Mach  numbers.  If  the  body  has  a  pointed  nose,  or  if  there 
are  no  canards,  or  if  the  tails  arc  absent,  the  above  time  and  cost  can  <.<e  reduced 
considerably.  For  example,  consider  the  pointed  body-tail  configurations  of  Figures 
29  and  30.  These  take  less  than  two  minutes  of  execution  time  for  ten  Mach 
numbers  and  cost  about  S25.  it  is  thus  believed  that  the  present  method  is  very 
cost  effective  compared  to  experimental  method,  in  that  reasonably  accurate  results 
for  forces  and  moments  can  be  obtained  at  a  relatively  small  cost. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


1.  A  general  method  lias  been  developed  consisting  of  several  theoretical  and 
empirical  procedures  to  calculate  lift,  drag,  and  center  of  pressure  on 
wing-body-tail  configurations  from  Mach  number  zero  to  about  three  and  for 
angles  of  attack  to  about  fifteen  degrees. 

2.  Comparison  of  this  method  with  experiment  for  several  configurations  indicates 
that  accuracies  of  ±10%  can  be  obtained  for  force  coefficients  of  most 
configurations.  This  is  at  a  cost  of  S75  or  less  for  ten  Mach  numbers  in  the 
range  0  <  Mw  <3. 

3.  First-order  theory  can  be  used  in  conjunction  with  modified  Newtonian  theory, 
to  calculate  wave  drag  on  wings  with  blunt  leading  edges.  However,  at  low 
supersonic  Mach  numbers  (Mm  <  1.5)  relatively  large  discontinuities  in  pressure 
coefficient  exist  at  the  match  point. 

4.  Wave  drag  of  wings  with  variable  airfoil  section  and  thickness  to  chord  ratio 
along  the  span  can  still  be  calculated  by  conical  flow  theory  if  the  local  source 
strength  is  computed  at  each  point  on  the  wing  based  on  the  local  slope  m 
the  point. 

5.  Changes  in  body  base  pressure  due  to  the  presence  of  tail  surfa- „s  is  not  known 

sufficiently  for  all  airfoil  shapes.  It  is  recommended  that  a  systematic  wind 
tunnel  study  be  conducted  for  Mach  number,  0  <  <  3,  angle  of  attack, 

0  <  a  -<  20°,  and  various  airfoil  geometries  to  empirically  estimate  this  change. 

6.  Although  much  progress  lias  been  made  in  transonic  flow  theory,  no  simple 
method  exists  for  calculating  wing  lift  or  wave  drag  on  complicated  wing 
configurations.  It  is  thus  recommended  that  much  work  continue  in  this  area 
until  practical  tools  are  developed. 
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appendix  a 


GEOMETRY  FOR  BLUNT  LEADING  EDGE  WING 


Assuming  the  wing  has  a  cylindrical  leading  edge  when  viewed  in  a  plane 
normal  to  the  wing,  it  will  have  an  elliptical  cross-section  when  viewed  in  a  plane 
parallel  to  the  frccstream.  The  equation  of  an  ellipse  is: 


— ‘  —  =  l 

ai  4 


(A-l ) 


where  a,  and  a,  are  the  semi-major  and  semi-minor  axis,  respectively.  If  the  radius 
of  the  cylindrical  leading  edge  is  denoted  by  rLL;  then: 


a,  =  rL[i./cosA, 


(A-2) 


a2  rLE 


Differentiating  Equation  (1),  the  local  slope  at  any  point  along  the  elliptical  airfoil 
leading  edge  is: 


cos  Aj^ 
z 


x  cos  A . 


<A-3) 


Defining  5  as  the  angle  which  the  tangent  to  the  surface  at  the  point  x.y.z  makes 
with  the  x  axis,  then 


6(x,y,z)  -  tan  1 


(A4) 


Here  the  leading  edge  radius,  rL ,  is  allowed  to  vary  in  the  spanwise  direction. 

Now  if  it  is  desired  to  determine  the  coordinates  of  a  given  point  on  the 
ellipse  given  a  value  of  6,  then  Equations  (A-l)  and  (A-2)  can  be  solved 
simultaneously  to  obtain: 


A-l 


X 


(A-5) 


=  rt  r, 

cos  A 


=  rLE(y)  cos  A, 
tan  5 


^  xcosA, 


rLF.(V) 


(A-6) 


where 


F  = 


( cos  A,\  2 
^tan6  / 

i  +  2 

tan  6 


APPENDIX  B 


NEWTONIAN  WAVE  DRAG  COEFFICIENT  OF  BLUNT  LEADING  EDGE  WING 


The  wing  plan  form  and  coordinate  systems  are  shown  in  the  above  sketch. 
Consider  a  cross-section  perpendicular  to  the  leading  edge  or  along  the  x'  axis.  The 
leading  edge  will  be  cylindrical  in  shape  as  seen  below. 

Z 


B-l 


The  elemental  force  acting  on  the  surface  area  of  the  wing  leading  edge  is: 


dCx<  = 


where 


P  cos  0 

1 


dA 


(B-l ) 


V-  S  f 

2  ^ °°  «  ref 


dA  =  rLE(y  )d0dy 
Hence,  for  one  fin  we  have 


e''  *  s 


f  n 


0\\  /•!>  1(2  cos  A  I  ) 


G 


p_  p 


cos  Or,  ,;(y')dOdy' 


-p  V- 

•y  ~  oo  oo 


(B-2) 


rc 


b/(2  cos  a  | ) 


cos  0rLE(y')d0dy' 


But  the  last  integral  of  Equation  (B-2)  is  zero  when  considered  over  the  entire 
airfoil  surface  so  Equation  (B-2)  may  be  written: 


• =  i  I 


o  b/(2  cos  A  j . 


Cp  cos  0rL).  (y')dOdy'  (B-3) 


o  o 

For  modified  Newtonian  theory. 


Cp  =  Cp  cos2  0  (B-4) 

1  1  o 

where  0  is  the  angle  between  a  normal  to  the  surface  and  the  freestream  velocity 
vector. 


B-2 


To  find  0  define  t?  as  the  unit  vector  normal  to  the  surface  so: 


A 

n  - 


(B-5) 


Here  7jy..  r)?,  are  the  direction  cosines  in  the  £x,,  $7.  directions, 

respectively  and  arc  =  cos  G .  7?y,  =  =  sin  0.  Now  according  to  the  dot 

product  of  two  vectors.  A.*  13  =  |A|  |B|cos0,  one  may  write: 


cos0  = 


V  ‘ft 
»  * 


IV, 


(13-6) 


But, 


V 

oo 


U  C,  ,  T  V 


e, ,  +  w'  (■- 

y  »  z 


(13-7) 


where 


tlL  =  V,  cos  .A, ,  /  -  V,  sin  A, ,  and  =  0. 

Substituting  (B*5)  and  (B-7)  into  (B-6)  and  performing  the  dot  product 
there  is  obtained: 


cos  0  ~  cos  Aj  cos  0 


Cp  -  fp  cos*  A,  cos2  0 


Thus, 


(13-K) 


and  (B-3)  becomes: 


APPENDIX  C 


TRAILING  EDGE  SEPARATION  DRAG 


The  incremen'al  force  on  the  elemental  area  of  the  wing  trailing  edge  shown  in 
the  sketch  below  is: 


In  analogy  to  the  work  of  Reference  1,  an  empirical  expression  which  takes  into 
account  the  slope  of  the  wing  trailing  edge  is 


Again,  assuming  a  linear  variation  of  the  trailing  edge  diameter  from  (dR)T,  at  the 
root  to  (d1)Tl  at  the  tip,  the  diameter  at  any  y  location  is: 


dTl  (y)  -  idR  )T  j 


(dT  )TI. 


!r ’t  i] 


where 

y  =  2y/h. 


(-1 


Then  I  quation  (C-!)  becomes: 


cA  =  ^Xbf  dTl;(p> 


n  2S 


WTt(y)  dy  = 


bCn  ri 

■tl 


dTK(y)dy 


Substituting  (C-3)  into  (C-4)  and  integrating  yields 


hCpu  r 
=  35-  [IR> 


k  ’  r  !■:  +  ^  • 


or  for  four  fins: 


C,  =  — 
'  B  S 


pii  r  i  -bop  _ 

7  =  y~  [(RrVe+<Rt)J 

J  ret  u  J 


APPENDIX  D 
GLOSSARY 


n 


GLOSSARY 


m 

b 


c 


Aspect  ratio 

Wing,  span  (does  not  include  body  radius) 

Chord  length  at  any  point  along  span 
Axial  force  coefficient  =  Cn 
Axial  skin-friction  drag  coefficient 
Trailing  edge  separation  drag  coefficient 
Zero  lift  drag  coefficient 
Wave  drag  coefficient  of  wing 
Cross-flow  drag  coefficient  of  wing 
Spanwisc  wave  drag  coefficient 
Laminar  ski.. -friction  coefficient 
Turbulent  skin-friction  coefficient 

Pitching  moment  coefficient  measured  about  nose  tip  (positive  nose 

Spanwise  pitching  moment  of  wing 

Pitching  moment  coefficient  derivative  (dCM/da) 

Normal  force  coefficient 

Normal  forej  coefficient  on  tail  of  body  shed  vortex. 

Normal  force  coefficient  on  tail  due  to  canard  shed  vortex. 

Normal  force  co.ffic.ent  derivative  (dC’N/da) 


D-i 


GLOSSARY  (Continued) 


Cn  Spanwise  normal  force  of  wi  ig 

Cp  Pressure  coefficient 

C„  Base  pressure  coefficient 

1  u 

C„  Stagnation  pressure  coefficient 

po 

cri  Chord  length  at  root  of  the  leading  edge  portion  of  a  modified  double 

wedge  airfoil  section 

cr2  Chord  length  at  root  of  the  trailing  edge  portion  of  a  modified  double 

wedge  airfoil  section 

ACp  Lifting  pressure  on  wing 

d  Body  diameter 

E,F,K  Elliptic  integrals 

f.  Canard  vortex  semispan  at  canard  trailing  edge 

h  Height  of  canard  shed  vortex  above  body  axis  at  tail  center  of  pressure 

i  Tail  interference  factor 

KB  Ratio  of  lift  of  body  to  that  of  wing 

KB(w)  Ratio  of  lift  of  body  in  presence  of  wing  to  that  of  the  wing  alone 

Kc  Ratio  of  lift  oi  canard-body-tail  combination  tc  that  of  wing  alone 

KVV([S)  Ratio  of  lift  of  wing  in  presence  of  bodv  to  that  of  wing  alone 

k,mv)  Ratio  of  lift  on  body  in  presence  of  wing,  due  to  a  wing  deflection  6,  to 

that  of  wing  alone 

k(  kj  =  tanAj 
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GLOSSARY  (Continued) 


kW(i>)  ^at‘°  oi  w'nS  in  presence  of  body,  due  to  a  wing  deflection  5,  to 

that  of  wing  alone 

Lc  Lift  of  canard-body-tail  combination 
M  Mach  number 

m  m  =  (J  cot  A, 

Rn  Reynold’s  number 

r  Radius  of  body  (variable) 

rLJ;  Radius  of  wing  leading  edge  in  a  plane  normal  to  leading  edge 

r,.R  Radius  of  wing  trading  edge  in  a  plane  normal  to  trading  edge 

Srcf  Reference  area  (maximum  body  cross-sectional  area  urn  css  wing  alone  is 

considered  in  which  case  the  wing  planform  area  is  used) 

S  Wing  planform  area 

s  Wing  semispan  plus  bod.v  radius  at  wing  root  chord 

Tw  Wall  temperature 

t  Wing  thickness  (variable) 

u,v.w  Perturbation  velocities  in  the  x,y,z  directions,  respectively 

V  Total  velocity,  V2  =  (U^  +u)2  +  v2  +w2 

x,y,z  Rectangular  coordinate  system  with  x  at  nose  tip,  y  out  right  wing,  and  z 
positive  up.  If  no  body  is  present,  x  begins  at  wing  root  chotd. 

o  Angle  of  attack 

0  VSTT 


D-3 


GLOSSARY  (Continued) 


r  Circulation  at  wing-body  juncture  of  combination 

7  Ratio  of  specific  heats  (7  =  1.4  for  air) 

5  Wing  deflection  angle 

i)  r,  =  k.Jp  =  I/m 

0  Angle  between  a  tangent  to  the  local  wing  surface  and  freestream  direction 

0]  Angle  on  blunt  leading  edge  where  Newtonian  theory  stops  and 

perturbation  theory  begins  (match  point) 

Aj  Sweepback  angle  of  a  wing  generator  (i  =  1,2, 3,4)  with  i  =  1  the  wing 
leading  edge  and  i  =  4  the  wing  trailing  edge 

X  Ratio  of  tip  cnord  to  root  chord  (etA.f) 

ja  Mach  angle,  /x  =  sin" 1  ( 1 /Mw  > 

v  Coefficient  of  kinematic  viscosity 

a  0  =  kj.v/x 


Velocity  potential 


X 

Wedge  half  angle  (measured 

parallel 

to  freestream)  of 

Subscripts 

a 

Conditions  at  transonic  Mach 

number 

Ma  =  Mfb  +  0.07 

b 

Conditions  at  transonic  Mach 

number 

Mb  =  Mfb  +0.14 

cp 

Center  of  pressure 

fb 

Force  break  Mach  number 

IA 

Infinite  afteibody 
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Subscripts 
C  Laminar 

NA  No  afterbody 

r  Root  chord 

SA  Short  afterbody 

T  Turbulent 

t  Tip  chord 

VIS  Viscous  component 

I  Wing  for  which  slender  body  values  of  interference  lift  arc  known. 

II  Wing  with  sweep  for  which  representative  values  of  interference  lift  are 
desired 

■»  Frcestrcam  conditions 
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